% procedure that calculates 1st derivative of time-varying elastance

dEdt_LV=gradient(E_t_LV)./gradient(time);




